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I. INTRODUCTION 


A. BACKGROUND 

The importance of environmental effects on communication links and radar systems 
has been recognized for the last few decades. Low altitude electromagnetic wave propagation 
over the ocean is strongly dependent on atmospheric conditions and sea state. It is essential 
to accurately characterize atmospheric turbulence in order to predict the propagation of 
electromagnetic signals. For a signal skimming the ocean surface, an understanding of ocean 
surface wave behavior should be given equal emphasis in electromagnetic wave propagation 
modeling. Specifically, it is necessary to relate the variation of the index of refraction to 
atmospheric turbulence and sea state. 


The refractive index at optical frequencies can be approximated by 


n-1=77.6(1+7.52x10 °A*)(P/T)10°°, (1) 


where P is the atmospheric pressure in millibars, 7 is the temperature in Kelvin and A is the 
wavelength of light in um. To determine the variation in the index of refraction, it 1s sufficient 
to know the variation of temperature in the atmospheric surface layer. In the surface layer, 
the vertical profiles of wind speed, u, temperature, 7 and humidity, g, all have logarithmic 


forms & =(u, 7, q), given by [Ref. 1] 


E=E, +E, [In(2/zo.) We (O)VK , (2) 


where €. is the flux scaling parameter, €, is the value at the surface, Zoe 18 the roughness 
length, « is the von Karman constant, ¢ = z/L, L_ the Monin Obukov length and wy, the 
empirically derived stability correction functions. 

Equation (2) shows that the model results in large refractive index gradients in regions 
close to the sea surface. An experiment conducted by the Naval Research Laboratory in June 
1994 [Ref. 2] showed that the result of the large index gradient caused rays from the 
transmitter to be bent more than that required to establish optical continuity between the 
transmitter and the receiver. 

Unlike the atmospheric boundary layer above land, the boundary layer above the sea 
is continuously interacting with the ocean aarti waves. The temperature profile of Equation 
(2) derived for an immovable surface may not be valid in the region of a few meters above the 
sea surface. Experimental data on signal heights that commensurate with ocean wave heights 
demonstrate deviations from the logarithmic distribution given by Equation (2) [Ref. 3]. 
B. OBJECTIVES OF THE RESEARCH 

This research evaluates the effects of ocean surface waves on optical fluctuations and 
relates the optical refraction to the near-surface gradient of temperature and ocean surface 
waves. Measurements from an optical link are used to infer the temperature profile and hence 


the refractive index profile. These are used to assess the effect of the ocean waves. Based on 


the experimental data, a proposed temperature profile is used in conjunction with a ray tracing 
program to validate optical continuity, | 
C. ANALYSIS APPROACH 

The research will first examine the small-scale atmospheric fluctuations due to the 
presence of the ocean surface waves on the optical signal fluctuations. Computed ocean 
surface spectral density and air temperature fluctuation spectral density are used to establish 
the influence of the ocean surface waves. Having established the influence of the ocean waves, 
modification to the near-surface mean temperature gradients is proposed. This is done by 
incorporating a random scaling length to account for the random motion of the ocean waves. 
With the proposed temperature profile, a ray tracing code was used to establish optical 
continuity between the transmitter and the receiver. 

This report comprises two volumes. Volume I contains seven chapters covering the 
following details: The experimental set-up and instrumented flux buoy are described in 
Chapter II. Chapter III presents the ray tracing code used in conjunction with the proposed 
temperature profile to validate optical continuity between the transmitter and the receiver. 
The analysis of ocean surface wave spectral density is presented in Chapter IV and the bulk 
aerodynamic method used to analyze the meteorological data is covered in Chapter V. The 
results of data analysis are presented in Chapter VI. Chapter VII summarizes the findings and 


recommendations of this research. The complete set of data plots is contained in Volume II. 





IL. EXPERIMENTAL SET-UP — 


A. OPTICAL LINK 

In the experiment, an optical link at 0.83 im was set up from Point Pinos at 
Lat./Long. (36° 38 10.931 N, 121°55 45.444" W) to Marina State Beach at Lat./Long. (36° 
42 01.365 N, 121° 48 33.651" W) . The propagation path was approximately 12.8 km and 
the bearing of the transmitter from the receiver was 237°. A sketch of the optical path is 
shown in Figure 1. 

The transmitter consists of a 0.83 um GaAs laser source with an average power of 
35mW and intensity modulated at 4.2 kHz. It has a beam aperture of 50 mm and a beam 
divergence half angle of 1.3 mrad. The receiver has a 50 mm aperture, a laser line filter with 
a bandwidth of 40 nm centered at 0.83 um followed by an Oriel Instrument Inc. battery 
operated silicon photodetector and pre-amplifier assembly. Block diagrams of the transmitter 
and receiver are shown in Figures 2 and 3 respectively. During the experiment, the heights 
of the transmitter and receiver were varied and the received signal recorded on tape. 

The recorded signal was demodulated with a Stanford Research Systems lock-in 
amplifier with sensitivity set at 50 mV, time constant of 1 ms and internal oscillator set at 4.2 
kHz. The output of the lock-in amplifier was connected to a Labview National Instruments 
12 bit analog to digital card plugged into a Macintosh computer. The analog signal was 
sampled at 400 Hz. A total of 41 signal samples covering the period from 30 May 96 to 5 Jun 


96 were taken. 





Marina (36° 42 01.365" N, 121° 48 33.651" 


Buoy (36° 40’ N, 121° 52 W) 


Point Pinos (36° 38 10.931" N, 121° 55' 45.444" W) 


Figure 1. Optical Path across the Monterey Bay 
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B. FLUX BUOY 

The flux buoy was deployed at approximately the mid point of the optical link, at 
Lat./Long. (36° 40 N, 121°52 W) ie the ocean depth is about 69.5 m. The flux buoy was 
deployed to record meteorological and wave data. The flux buoy system contains two 
separate sensor suites, viz., a mean meteorology package and a high frequency sampling ‘flux’ 
package. Two sets of data were collected; the high frequency data (flux) of air temperature 
and wind velocity fluctuations sampled at 21 Hz, wave motion data sampled at 10.5 Hz and 
low frequency data (mean meteorology) of wind speed, wind directions, humidity, pressure, 
air and sea temperature sampled at one minute interval. A picture of the buoy is shown in 
Figure 4. 

The high frequency air temperature and velocity were measured by the Solent 
ultrasonic anemometer at 5.09m above mean tide level. The measurement technique is based 
on sound pulses transmitted back and forth between transducers. The time of flight in each 
direction, say ¢, and t, are measured. If c is the speed of sound, x the distance between the 
transducers, and there is an air flow v along the line of the transducers, the following 


relationships are readily derived: 


Sy eiareriae (3) 


i = (4) 


(c-v) © 


1 
v=05x i | (5) 


and 


C= 0) x¥ 





fo 
Ze (6) 


I 2 


By arranging three pairs of transducers in different orientations, the direction and magnitude 
of the incident air flow may be unambiguously derived. The air temperature is calculated from 


the speed of sound as follows [Ref. 4]: 


c* = yR7(1+0.51q) , (7) 


where yR is 403 m? s* K" for air, 7 is the air temperature in degrees Kelvin and gq is the 
specific humidity in kg of water vapor per kg of air. 
The ocean surface wave data was measured by the Systron Donner Motion Pak. The 


Motion Pak consists of three accelerometers measuring the accelerations of the x, y and z 
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directions (sway, surge and heave) and three rate gyros measuring the angular rates, 8, o and 
W (pitch, roll and yaw). The data was sampled at 10.5 Hz. The Motion Pak uses a night- 
handed coordinate system for the angular rates about the x, y, and z axis. It was mounted so 
that the x-axis points towards the port side of the buoy, the y-axis points towards the bow and 
the z-axis downward. The z-axis voltage was adjusted so that when the buoy is at rest on a 
flat surface, the z-axis linear acceleration voltage output read zero voltage. 

Three wave staffs were deployed to measure the high frequency ocean surface waves. 
The three wave staffs were arranged such that each staff Res located at the vertex of an 
equilateral tnangle. The arrangement allows for the computation of the wave directional 
spectra. Unfortunately, wave staff number 2 was not functioning properly during the 
experiment. 

The buoy’s heading relative to magnetic North can be determined from the three 
Syntron compass outputs. The magnetic deviation from the true North at the buoy was 15° 
15 East. The buoy’s heading relative to the true North can be determined by subtracting the 
deviation from the buoy magnetic heading. 

The low frequency wind speed was measured by the R.M. Young Wind Monitor 
located at 3.80 m above mean tide level. The wind speed sensor is a helicoid shape propeller 
molded of polypropylene plastic. Threshold sensitivity of the propeller is 1.0 m/s. Rotation 
of the propeller produces a linear DC voltage signal directly proportional to the wind speed. 

The Everest IR sensor was used to measure the sea surface temperature. The IR 


radiation from the sea surface 1s collected by an optical lens system and focussed on an IR 


1] 


detector. The detector in turn converts the IR signal to an electrical signal proportional to the 
sea surface temperature. 

The low frequency temperature and humidity were measured by the Rhotronics 
relative humidity and temperature sensor. The humidity sensor is a small hygroscopic 
capacitor that changes its value as a function of both water vapor pressure and temperature. 
The temperature sensor uses a platinum 100 ohm resistance temperature device (RTD) to 
measure air temperature. 

The flux data were stored in an Onset Model 6 Toshiba computer with 500 MB of 
disk storage on board the buoy. Each word is a two-byte data stored as ‘big endian’ 
hexadecimals (most significant byte first). The mean meteorology data were stored in another 
Onset Model 6 computer with 200 MB of disk storage. The data were stored in ASCII 
format. The bulk aerodynamical method was used to — the mean meteorological sata 
Cc; ENVIRONMENTAL INFLUENCE AND LIMITATIONS 

The Monterey Bay has deep water with ocean-like swells. The optical link from Point 
Pinos to Marina State Beach crosses Monterey Bay where there is relatively very little 
shoreline influence. The buoy was located at the mid point of the link in the deepest water and 
‘farthest from land with open-ocean wind and waves. 

The limitation in the deployment of the buoy at a single location leads to the 
assumption of a spherically stratified atmosphere with a lateral homogeneous refractivity 
along the path of propagation. This assumption is reasonable only in a well-mixed 


environment over the ocean. 


In the experiment, the spatial average optical measurement is compared with the 
temporal average meteorological data. In a situation where there is lateral homogeneity along 


the mean wind direction, the two averages could be considered to be equivalent [Ref. 1]. 
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Figure 4. Flux Buoy 
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If. RAY TRACING 


The ray tracing code was used in conjunction with the standard logarithmic 
temperature profile and the proposed modified temperature profile to trace rays emanating 
from the transmitter to the receiver. Optical continuity is said to be achieved if there exists 
at least one ray above the ocean wave height connecting the transmitter to the receiver. 

The ray tracing technique used here is based on the methodology developed by Lehn 
[Ref. 5]. The model assumes that the atmosphere consists of a number of concentric spherical 
layers where the temperature is assumed to vary linearly within each layer. In this case, the 
curvature of a nearly honzontal ray 1s approximately constant within each layer. Replacement 
of circular arc by parabolic arc leads to a simpler quadratic expression for the ray trace. 

For a ray traversing a stratified atmosphere above a flat earth, Snell’s Law dictates 


that 


nsinO = A, (8) 


where 77 is the refractive index of the air at any point on the ray, 9 is the angle between the 
ray and the vertical, and A is a constant for the ray. 


In rectangular coordinates, the curvature, C, is given by — 


ts) 


(9) 


where x is the horizontal distance and z is the elevation above the earth’s surface. The slope 


of the ray is given by 
dz _ cosO 


dx sin® 
(n 2° A aye 
= — 





(10) 


and the second derivative is given by 


dz d ao 
2 dz| dx \dx 
a (11) 


With Equations (10) and (11), the ray curvature becomes 


(12) 
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The sign in Equation (12) has been chosen so that a positive curvature corresponds to a ray 
that is concave towards the earth. At optical frequencies, the refractive index of air can be 


expressed as 


Me | ee 
if 


; (13) 


where p is the pressure, 7 is the absolute temperature in degrees Kelvin, € = 226x10* m’ kg” 
and B = 3.48x10° J’ kg K. Note that Equation (13) is equivalent to Equation (1) when the 
dependence on wavelength is set at A = 0.83 ym and the pressure expressed in N/m’. The 


derivative of 7 is given by 


(14) 


where the relation 8 = pp/T has been used, p is the density of the air and g is the acceleration 


due to gravity. Substituting Equation (14) into the expression for the ray curvature leads to 


ae sind eBp Pepe 2 


(15) 





pe. a 
(1+€p) Ip _ | 
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In ray tracing, the circular arcs that represent the earth’s surface, the ray and the layer 
boundaries are all approximated by parabolic arcs. The coordinate system is shown in Figure 


5. The surface of the earth is given by 


2 eae (16) 


where R is the radius of the earth. Similarly a layer located at a height 4 above the earth’s 


surface is given by 


2 
x 
Zz +h 


2(R+h) 





(17) 
2 
ey, 
2R 


R 


An optical ray with a radius curvature, r = //C, starting at x = 0, z = z, and making an angle 


of , with the horizontal is given by 


2 
<s : 
z= -— + xtanod + Ze. 18 

2r c : (18) 
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Within the 7” layer, the ray path is given by Equation (18) with a radius of curvature 
r, . Ray tracing is performed by determining the intersection of the ray with the layer 
boundaries. The 7 layer lies between elevations /., and h,. A ray can enter a layer either from 
above or below and could exit the layer either from top or bottom of the layer. Figure 6 
shows a ray entering the i” layer from lower boundary. It will proceed to P ifr, > R or r, < 
0 and to Q if r, is positive and sufficiently small. The coordinates of P are found by solving 
Equations (17) and (18) with A = 4, in Equation (17) and z, = A, in Equation (18). The x- 


coordinate of P is 


Rr 


i 


rik 


Up = 





(19) 








-tand, tan’, a = mal = 


mK 


If the discriminant is negative, the intersection P does not exist. Then Q has to be determined. 


This is done by letting # = h,, in Equation (17). The x-coordinate of Q is 


2Rr,tang, 


0) ee 
Q Ror (20) 


Sumilarly, for a downward-heading ray, the x-coordinate of P at layer boundary h,, is given 


i 


by 


1/2 


— 24,-h,.0,-P) 


OR (21) 





Up = — -tand, -| tan’, 
PAR 





Again, a negative discriminant implies that P does not exist and the ray bends upwards and 
returns to layer boundary h, at Q. The x-coordinate of Q 1s determined from Equation (20). 

New coordinates are used each time the ray crosses a layer boundary. Figure 7 shows 
the case of an upward-heading ray. As the curvature is small, the x-coordinate of the new 
origin is taken to be wp. The layer boundaries can still be described by Equation (17). The 


elevation angle of the ray is given by 


od, = tan'm, - tan'm, , (22) 


where m, and m, are the ray and layer slopes respectively at P. The slopes are found by 


differentiating Equations (17) and (18) giving 


R 


tanh, = m, - m, 
(23) 
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Figure 5. Ray Tracing Coordinate System [From Ref. 5] 





Figure 6. Upward-heading Ray in Layer i [From Ref. 5] 
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Figure 7. Change of Coordinates when the Ray enters a New Layer [From Ref. 5] 
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IV. OCEAN WAVE SPECTRAL ANALYSIS 


A. MEASUREMENT APPROACH 

In order to establish the possible effects of ocean surface waves on the optical link, 
it is necessary to examine the ocean surface wave spectra and wave heights at the time when 
each sample of optical signal was measured. The wave information was collected by the 
motion pack and wave staff. The motion pack consists of three accelerometers measuring the 
accelerations of the x, y and z directions (sway, surge and heave) and three rate gyros 
measuring the angular rates, 8, @ and w (pitch, roll and yaw). The data was sampled at 10.5 
Hz. 

The true surface displacement is obtained by summing the following three 
components, viz., the wave height as measured by the wave staff, z,, the displacement of the 
wave staff due to the translational motion of the buoy, z, , and the displacement of the wave 
staff due to the rotational motion between the wave staff and the motion pack [Ref. 6]. 

As the wave staff and the linear accelerometers were mounted on the moving 
reference of the buoy, it is necessary to compute the surface displacement in an earth- 
referenced system. The conversion of the buoy-referenced system to the earth-referenced 


system can be achieved via a transformation matrix, T,; [Ref. 7] given by 
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cosOcos sindsinOcospy-cospdsiny  cosdsinOcosy + sind sinyy 
= |cosOsiny  sindsinOsinb+cosdcosy cosdsinOsiny -sindcosy|. (24) 


-sin8 sind cosO cos®cosO 


T 


BE 


If a, represents the translational accelerations measured by the strap-down 


accelerometers on the buoy, the earth-referenced accelerations, a, are given by 


a,=T1,,4,+8, (25) 


where g = (0, 0, -g) is the acceleration due to gravity. The earth-referenced vertical 


displacement due to the translational motion of the buoy is given by 


ee [ [ [a, . (-sinO, cosO sing, cosO cosh) - g]dtdt . (26) 


Similarly, the earth-referenced vertical displacement measured by the wave staff 1s given by 


Zz = ZmpC0SO cosh . (27) 
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Finally the displacement due to the rotational motion is determined by integrating the 
rotational velocity vector, Q x T,; L, where Q is the angular velocity of the buoy measured 
by the rate gyro and L, is the distance between the accelerometer and the water surface at the 


wave staff. The angular velocity, Q is given by 


-O sin + cos cosy 
Q = | Ocosp +cosOsiny | - (28) 
y- sind 


The vertical displacement due to the rotational motion is given by 


— [Le2(-8 simp +cosOcos)-L,,(Ocosy +cosO sin) Jat , (29) 


where 


BE Lp (30) 


and the subscripts £/ and £2 represent the component of the vector. 


The true ocean surface displacement is given by 


2a 


22 ee (31) 


Once the ocean surface displacement is determined, the mms ocean surface displacement, o 


can be computed and the significant wave height, H, is given by 


H_ = 40. (32) 


B. ANALYSIS OF OCEAN SURFACE DISPLACEMENT 

The analysis of ocean surface displacement or elevation was carried in the frequency 
domain using Fourier Transform. The rate gyros do not have good low frequency response 
and therefore it 1s necessary to use con Memehtary filtering to patch the low frequency tilt 


reference of the rate gyros with the linear accelerometers as follows 


0 = 6.,-sin '(@,/g),> » (33) 


b = $+ sin ‘(a,/g), > (34) 


and 


y = W,p+sin ‘(com),, , (35) 


where a,, and a, represent the observed linear acceleration in the x and y direction, the 
subscripts HP and LP represent the high pass and low pass filtered output and com is the 
compass output. The first term on the nght hand side of Equations (33) to (35) represents the 
high pass angular motions from the rate gyros which provides the wave-induced angular 
motions and the second term represents the low frequency tilt reference from the 
accelerometers. The effect is an all-pass filter which removes the unwanted drift in the rate 
gyros while retaining the low frequency tilt reference. 

To compute z ,, and z,,, the integration and double integration are done in the 
frequency domain by dividing the Fourier Transform of z ,, and z ,- with (iw) and (iw)? 
respectively. The division by w introduces a singularity at w = 0. However, there is no ocean 
wave with frequency lower than 0.04 Hz, 1.e., 25 second waves, since there is not enough 
fetch in the ocean to sustain the waves without breaking. The waves in the Atlantic Ocean and 
Pacific Ocean rarely exceed 20 seconds and 25 seconds respectively [Ref. 8]. The singularity 
problem 1s alleviated by performing a high-pass filtering of the wave spectrum at 0.04 Hz. The 
power spectral density of the ocean surface waves and significant wave height were computed 


for the duration of time when the optic signal samples were being measured. 
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V. BULK AERODYNAMICAL METHOD 


The bulk aerodynamical formulations for the fluxes of momentum, sensible heat and 


latent heat across the air-sea interface and vertical profile estimation are given by 


<= oC,(z) [u(z)-u J? : (36) 
O, = —pc,C,(z) [u(z)-u,] [7(2)-T,] (37) 

and 
OQ, = -pL,C,(z)[uZ)-u,][9@)-9,] . (38) 


where p is the air density, c, is the specific heat of air, L, is the latent heat of vaporization, Cp, 
C,, and C, are the height dependent bulk transfer coefficients for momentum, sensible heat 
ad latent heat respectively, u, 7, g are the wind speed, potential temperature and specific 
humidity respectively, z is the height above the sea surface and the subscnpt s denotes the 
respective values at the sea surface. 


In the atmospheric surface these fluxes can be expressed according to Monin-Obukov 
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Similanty (MOS) theory scaling parameters as 


T= pul, (39) 
Or a “pc,u, i (40) 

and 
o = -pLu.q, > (41) 


where u., 7’. and g. are the MOS surface layer scaling parameters for momentum, temperature 
and humidity respectively. By combining Equations (36) to (41), these scaling parameters can 


be defined according to the bulk transfer coefficients as 


u, = one [u(z)-u J , (42) 


fea ree ANZ) (43) 
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and 


qe= CPCs. le@)-4.)'. (44) 


The bulk transfer coefficients can be defined as 


2 


tp Se eet 
” [In(z/z,,) -W,(2/L)P (49) 
Colt) a (46) 
—— [In@/z,,,) - b(2/L) ][In(z/z,,) - wA2/L)] 
and 
Cz) = z (47) 


[In(z/2,,) - W(2/L) [in(z/z,,) - p(2/L)] 


where xk is the von Karman constant, and ,, W, and , are the integrated forms of the 


dimensionless stability correction functions for wind speed, temperature and specific humidity 
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respectively, given by 


_J  2in{(1+x)/2] + In[(1 +x7)/2] -2tan“'(x)+7/2, z/L<0 
wycety=| alt i 
2In[(1 +x)/2 /L<O | 
VoCelL)=w,Cell)=} MIO O7L ESO (49) 


x = (1-20z/L)"4 (50) 


and Zy, , Z9r and z,, are the roughness lengths for momentum, temperature and humidity 


respectively. The Obukov length scale, Z is defined as 


pe Tu, 
~ [ke(T. +0.617q,)] (S1) 


where g is the acceleration due to gravity and 7’, is the virtual temperature introduced to take 


into account the effects of humidity: 


T= 701 +0.61q) . (52) 


To determine the vertical profile of wind speed, temperature and specific humidity, 
iterative computation of u., Ts, Qs, Zp, Zr, 29g and L are performed and the vertical profile 


calculated as 


u=u,+ = [In(z/z,,)-W,@/L)] , (53) 
r, 
P= 7, +— [In (2lz9p) - Up (Z/E)] (34) 
and 
q=4,+~[hn (z/z,,)-W,(e/L)] (55) 
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The temperature profile estimated in this manner was used as input to the ray tracing 
program to determine if optical continuity between the transmitter and the receiver can be 


established. 
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VI. RESULTS OF DATA ANALYSIS 


A. OCEAN SURFACE WAVE SPECTRA AND WAVE HEIGHTS 

The power spectral densities of the ocean surface waves were computed for the 
duration when the optical signals were being measured. A total of 41 samples were computed. 
A typical sample of the wave spectrum is shown in Figure 8. It can be seen from the spectrum 
that the ocean surface waves are dominated by swells with a peak frequency of about 0.1 Hz, 
i.e. 10 second waves. There are wind generated waves with a peak frequency at about 0.4 Hz, 
i.e., 2.5 second waves. The spectrum exhibits a high frequency roll-off as w~. This is the 
established behavior of the equilibnum range of wind-sea spectra [Ref. 9 ]. 

The rms wave height, 0, was computed from the power spectral density by the 
following relation | 


ot = 1 y IZ(/)|? (56) 
N? i=l 


where Z(/) is the Fast Fourier Transform (FFT) of the ocean surface displacement, z(t) and 


Nis the sample size. For the 41 samples, the rms wave height ranges from 0.18 m to 0.7 m. 
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Ocean surface elevation spectrum, 30 May 96, 1225 hr 
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Figure 8. Ocean Surface Elevation Spectrum. The straight line represents the -4 slope. 
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B. TEMPERATURE SPECTRAL DENSITY 

The aspect of atmospheric turbulence of relevance to optical propagation is the 
variation of the index of refraction or equivalently, the variation of the density of air. Air 
density depends on pressure and temperature. Pressure fluctuations are rapidly dissipated and 
their effects are negligible compared to temperature fluctuations. Thus, for optical 
propagation purposes, temperature fluctuations are used to characterize the optical effects 
of turbulence. 

In the experiment, temperature fluctuations sampled at 21 Hz were measured by the 
sonic anemometer at 5.09 meters above mean sea level. A typical plot of the power spectral 
density of temperature fluctuations is shown in Figure 9. The spectrum shows the classical 
-5/3 slope at high frequencies. This dynamic condition in the high frequency range agrees 
satisfactorily with theoretical predictions and eo above land [Ref. 10]. At lower 
frequencies, the spectrum shows evidence of ocean wave influence. Figure 10 shows the 
temperature spectral density and corresponding ocean surface spectral density of the wind 
waves. The spectral density of wind waves was obtained by performing high pass filtering on 
the swell-dominated wave spectrum, thus removing the swells from the spectrum. It can be 
seen that the temperature spectrum displays the characteristic spectral density peak at the 


peak of the corresponding wind wave spectrum indicating the influence of the ocean waves. 
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Temperature spectrum, 30 May 96, 1225 hr 
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Figure 9. Temperature Spectrum. The straight line represents the -5/3 slope. 
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Temperature spectrum, 30 May 96, 1225 hr 
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Figure 10. Temperature Spectrum and the Corresponding Wind Wave Spectrum. 
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C: OPTICAL SIGNAL FLUCTUATIONS 

In the analysis of the optical signals, the signal fluctuations are described by the log- 
normal variance although the variance of the signals can be used as a measure of fluctuations. 
The research does not attempt to relate the log-normal variance to the saturation 
phenomenon. The objective is to analyze how the effects of ocean surface waves affect the 
signal fluctuations per se. 

1. Log-normal Distributions 

To observe the signal fluctuations, a histogram plot of the normalized log intensity 
(normalized to the mean intensity) was performed for each signal sample. From the histogram 
plot, it was found that for signal samples where no apparent ocean wave occultation is 
observed, the intensity fluctuations exhibit a log-normal distribution. For signal samples where 
wave occultation is evident, deviation from the log-normal distribution results. Typical 
histogram plots of samples that exhibit log-normal distribution and ‘non log-normal 
distribution’ are shown in Figures 11 and 12 respectively. 

Log-normal distributed signals are known to exist in signals propagating through 
random inhomogeneous medium [Ref. 11]. The distribution characteristics are associated with 
atmospheric radio noise or atmospherics as it is sometimes called. It has been found 
empirically [ Ref. 12] that the distribution P(E>R), where £ is the instantaneous value of HF 
atmospheric noise and & is threshold voltage at the receiving point, approaches a Rayleigh 
distribution for small amplitudes. This is due to the interferences of many mutually 


overlapping small atmospherics with uniformly distributed phase. However for large 
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amplitudes, the distribution deviates considerably from the Rayleigh distribution and 
approaches a log-normal distribution. This is believed to be due to strong atmospherics that 
do not overlap in time. 

The theory of the distribution of peak value of atmospherics can be explained as 


follows. The peak value of atmospherics, E, can be represented by the following expression 





By = ; on -)> sd (57) 


where P, is the peak power of the atmospheric discharge, dis the propagation path, d, is a 
section of the propagation path, 6, is the attenuation constant of the j” section and K is a 


constant. By writing Equation (57) as 


E, =e” = exp Le + InP, - In))d, + InK] , (58) 
jas 


the sum of the attenuations of each section, A, is defined. It is assumed that the fluctuations 
of P, and d are negligible compared to the first term in the exponential of Equation (58). 


Thus, only the first term will thus be used to determine the distribution of A. Since most 
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atmospherics are propagated from long distances away, it is reasonable to assume that the 
attenuations 6 }, d, are independent and none of the 6 ), d,’s make a predominant contribution 
to the sum of the attenuations. Therefore A is normally distributed with probability density 


function given by 








l A-py 
p(A) = ——exp Sb (59) 
O,y2T 20, 
where i is the mean, O% is the variance. From Equation (58) we have 
i (60) 
and substituting Equation (59) into Equation (60), the density function of £, is 
(InE,-p)? 
eet eee (61) 
0 EV2m 204 


Equation (61) is the log-normal distribution which has been found experimentally to apply to 
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the propagation of optical signals through atmospheric turbulence in the absence of 


occultation of the optical signals by ocean surface waves. 
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Optical signal, 30 May 96, 1225 hr 
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Figure 11. Optical Signal Intensity and Distribution Plots and the Corresponding Ocean Waves. 


‘.... IN Ocean waves plot indicates height of optical ray path above the waves. 
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Optical signal, 30 May 96, 1859 hr 
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Figure 12. Optical Intensity and Distribution Plots and the Corresponding Ocean Waves Plot 


(with Occultation of Optical Signal by Ocean Waves). 
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2: Log-normal Variance and Mean 
From the above equations, the relationship between the log-normal mean intensity and 
log-normal variance can be derived. In general, an electromagnetic wave field E(7,2) can be 


represented by 


E(t) = 1(7)'? exp[ik-? - iwt] , (62) 


where is the temporal frequency of the electromagnetic field, & is the propagation vector 


and /(7) is the irradiance or intensity. In the presence of turbulence, Equation (62) becomes 


EG) = <I(F)>'? exp[ik-F - iwt + x(%1) + 07D), (63) 


where<i/(7)>!” is the average non fluctuating part of the irradiance, y(7,t) describes the 
fluctuating part and (7,1) is the phase variation. The random variables, y(7,t) and (7,1) 
have Gaussian probability distributions. To examine the irradiance fluctuations, the 


instantaneous irradiance can be written as 


[(r,t) = <I(r)>exp[2x(7,0] . (64) 
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It follows that /(7,¢) has a log-normal distribution. Taking the average on both sides of 


Equation (64), and using the relation valid for any real-valued Gaussian random variable, 


(225 a ext i = a0? : (65) 





Equation (64) becomes 


<I> = <Izexp(2<x>+20,”) . (66) 


In the presence of turbulence, </> should be a constant independent of the turbulence. It 


follows that 


CE Boe (67) 


A plot of log normal mean intensity against the log normal variance for those signal 
samples not occluded by the ocean surface waves in Figure 13 further verifies the log-normal 


distribution characteristics of the optical signals. 
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Figure 13. Log-normal Variance vs. Log-normal Mean. 


48 


D. EFFECTS OF OCEAN SURFACE WAVES ON OPTICAL SIGNAL 
FLUCTUATIONS 


Due to the random motion of the ocean waves, it is necessary to establish a way to 
measure the distance of the ray path above the ocean waves. The method adopted here is to 
calculate the normalized mean square distance, /’, (normalized to the variance of the ocean 


surface displacement) of the ray path above the ocean waves, given by 


h? = “; Dd, Gz (68) 


where 0” is the variance of the ocean surface displacement, z,, is the lowest height of the 
straight line ray path above mean sea level, z, is the instantaneous ocean surface displacement 
and the summation is performed over the duration of the optical signal transmission. A large 
mean square distance implies a lower probability of the ray path being blocked by the ocean 
surface waves. The normalized mean square distance is plotted against the log-normal 
variance of the signals. For sufficiently large mean square distance, 1.e., for the mean square 
distance greater than 35, the variance of the optical signal increases as the mean square 
distance decreases. Those signal samples with lower mean square distance also display the 
same trend but with increased variability, 1.e., the signal sample points are widely scattered. 


A plot of log-normal variance of the optical signal against /’ is shown in Figure 14. The signal 
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samples are numbered from 3 to 43 for easy identification of any signal sample. The list of 
optical signal samples is contained in Appendix A. 

‘From Figure 14, it can be seen that there is an apparent ‘peak’ in the signal 
fluctuations for signal samples 14 to 17. The higher signal fluctuations are due to the 
increased roughness of the ocean wind waves associated with the signal samples. The wave 
age, defined as the ratio of the ocean surface waves phase speed, C, to the wind speed, w is 
a good measure of roughness of the surface waves. Donelan, et. al. [Ref 13], demonstrated 
that younger waves traveling slower than the wind are rougher than mature waves. A plot of 
the optical signal log-normal variance against wave age for signal samples 3 to 5 and 14 to 
17 is shown in Figure 15. The plot indicates a decrease in signal fluctuations with increased 
wave age. However, signal samples 14 and 15 have higher signal fluctuations than samples 
3 and 4 although they are associated with wind seule waves of approximately the ae 
wave age. The higher fluctuations of samples 14 and 15 are due to the influence of the swell. 
The influence of the swell on sample 14 and 15 is manifested in the temperature fluctuations 
spectral densities where the swell frequencies modulate the temperature spectra. This effect 
is shown in Figures 16 and 17. The corresponding temperature spectra for samples 3 and 4 
do not display any prominent peak in the swell frequencies. The temperature spectra and the 
ocean surface elevation spectra for samples 3 and 4 are shown in Figures 18 and 19. The 
presence of swell caused a change in the spectral distribution of the various components of 
the high frequency wind waves. There is experimental evidence that this leads to an increased 


in roughness of the ocean surface [Ref. 14]. 
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Two alternative analyses were performed to verify the trend of the signal fluctuations. 
In the first method, it is assumed that the ocean surface elevation follows a Gaussian 
distribution, n(0,0). This is true for a sample with a large number of-points. Based on the 
height of the ray path above mean sea level, the probability of ocean surface waves protruding 


into the ray path can be computed from 


oo 


P(z>2,,) = [ = e 2120") gy (69) 


2 
2, 2TtoO 


The percentage of the ocean waves protruding into the ray path is found by multiplying 
Equation (69) by 100. The percentage values are plotted against the signal log-normal 
variance. It was found that the signal fluctuations exhibit the same trend. A plot of the 
percentage of ocean waves protruding into the optical ray path is shown in Figure 20. 

In the second method, the log-normal variance of the signal samples are plotted 
against the difference between the ray path height and one significant wave height above the 
wave crests. The same trend in signal fluctuations was observed. For signal samples that are 
less than one significant wave height above the ocean surface waves, the sample points are 
widely scattered as in the above two analyses. The plot is shown in Figure 21. The analyses 


led to the conclusion that the random motion of the ocean waves is affecting air turbulence 


ul 


immediately above the ocean surface and causing increased variability in signal fluctuations 


as the ray path gets closer to the ocean surface. 
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Figure 14. Log-normal Variance of Optical Signal vs. h’. 
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Figure 15. Optical Signal Log-normal Variance vs. Wind Wave Age. 
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Temperature spectrum, 31 May 96, 951 hr 
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Ocean surface elevation spectrum, 31 May 96, 951 hr 
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Figure 16. Optical Signal Sample No. 14, Temperature Spectrum and the 


corresponding Ocean Surface Wave Displacement Spectrum 
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Pelle tie spectrum, 31 May 96, 1115 hr 
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Ocean surface elevation spectrum, 31 May 96, 1115 hr 
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Figure 17. Optical Signal Sample No. 15, Temperature Spectrum and the 


Corresponding Ocean Surface Wave Displacement Spectrum. 
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Temperature spectrum, 30 May 96, 1225 hr 
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Ocean surface elevation nspectmum, 30 May 96, 1225 hr 
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Figure 18. Optical Signal Sample No. 3, Temperature Spectrum and the 


corresponding Ocean Surface Wave Displacement Spectrum 
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Temperature spectrum, 30 May 96, 1230 hr 
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Ocean surface elevation spectrum, 30 May 96, 1230 hr 
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Figure 19. Optical Signal Sample No. 4, Temperature Spectrum and the 


corresponding Ocean Surface Wave Displacement Spectrum 
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Log-normal variance vs. percentage of ocean waves protruding into ray path 
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Figure 20. Percentage of Ocean Waves Protruding into Optical Ray Path. 
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Figure 21. Log-normal Variance of Optical Signal vs. Ray Path 1 Significant Wave Height 


above Ocean Waves. 
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E. OPTICAL SIGNAL MODULATION BY OCEAN SURFACE WAVES 

The intensity spectrum for each optical signal sample was computed and compared 
with the corresponding ocean surface wave spectrum. The plot (Figure 22) shows that for 
optical signals that were partially occluded by the ocean surface waves, the optical signal 
spectrum displays peaks at the swell frequency of the ocean waves. The optical signal is said 


to be modulated by the ocean waves. 
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Figure 22. Optical Signal Spectrum and the Corresponding Ocean Wave Spectrum. 
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F, NEAR-SURFACE MEAN GRADIENTS 

In the above, only small-scale variations are considered whereby turbulence produces 
random fluctuations in the density of air and hence fluctuations in the refractive index. On 
larger scale vanations, the main concern is on phenomena that affect the average position of 
the optical beam as it traverses the atmospheric medium. In the experiment, the effects of the 
ocean waves on vertical temperature profile and hence the refractive index profile is 
investigated. 

There is evidence from the experiment that the air turbulence due to the ocean waves 
resulted in deviations of the vertical temperature profile from the loganthmic profile thereby 
affecting the refraction of the optical signals. To analyze the refraction of the optical signal, 
only samples of those above the ocean waves, i.e., those not occluded by ocean waves, were 
used. For each time sample, the corresponding transmitter and receiver heights above mean 
sea levels were calculated from GPS surveying data and tide measurements provided by a 
National Oceanographic and Atmosphenc Administration tide gauge. With the measured 
transmitter and receiver heights, and using the standard logarithmic temperature profile in the 
ray tracing program, it was observed that optical continuity between the transmitter and the 
receiver could not be established. An example of the ray tracing plot is shown in Figure 23. 
The mechanism in which the ocean waves affect the index of refraction needs to be examined. 
This research attempts to relate the effects of ocean waves on the temperature profile using 
the mixing length concept. 


The mixing length concept was introduced by Prandtl in 1925 to describe the 
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properties of turbulent transfer. In this model, it is assumed that there exists a length scale 
called the mrng length over which a turbulent eddy can travel and still maintain its identity. 
Details of the mixing length concept is contained in Appendix B. 

The effects of ocean surface waves on the gradients of the vertical temperature profile 


are discussed in the next section. 
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Ray tracing using standard logarithmic temperature profile, 30 May 96, 1225 hr 
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Figure 23. Optical Continuity cannot be established between Transmitter and Receiver. 
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G. EFFECTS OF OCEAN SURFACE WAVES ON THE INDEX OF 
REFRACTION 


The refractive index at optical frequencies can be approximated by 


n-1=77.6(1+7.52x107A*)(P/T)10° , (70) 


where P is the atmospheric pressure in millibars, 7 is the temperature in Kelvin and A is the 
wavelength of light in pm. In the atmospheric surface layer, pressure variation is small and 
can be neglected. Variation in the index of refraction is thus caused by the variation of 
temperature in the atmospheric surface layer. 

The vertical profiles of wind speed, u, temperature, 7 and humidity, g all have 


logarithmic forms € =(u, 7, q), given by 





(71) 





and the gradient profile given by 
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- 0 z) ; (72) 


where , is the flux scaling parameter, €. is the value at the surface, z), is the roughness length, 
K is the von Karman constant and @,(z/L)is the derivative version of the stability functions 
given in Appendix B. In Equation (72), the term xz is the turbulence mixing length. 

The above profiles are derived for turbulence above an immovable surface. A 
distinctive feature of the boundary layer adjacent to the ocean surface is the presence of a 
continuous random movement of the ocean surface. The presence of the random varying air- 
wave interface produces a constant variation in resistance of the ocean surface. The 
roughness elements that originate from the various components of the wave spectrum 
contribute in varying degree to the total wave dispersion. The air flow above the ocean waves 
induces a pressure field which transfers part of the energy of the air flow to the waves during 
wave development or conversely, the waves transfer energy to the air during swell 
attenuation. The surface waves, on the other hand, disturb the air flow above them and form 
air waves which interact with the gradient of the vertical profile. 

The wavy motion of the surface affects the profile deviation in two ways. Firstly, the 
short waves (high frequency component) contribute to the roughness of the water surface. 
The quantity z,, specifies the scale of the turbulent eddies that are generated by the surface 
roughness elements. Experimental evidence [Ref. 16] suggested that surface roughness should 


be proportional to wave steepness, H/A, where H is the wave height and A is the wavelength 
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of the dominant ocean wind waves, 1.e., 








Pie 
Pa! 73 
Ou A g ( ) 
As the ocean wave phase speed, C « [gA/(27)]'”, z,, can be expressed as 
au.) 
2. = (74) 
= a c 7 


™ 


where a is the constant of proportionality. Equation (74) implies z), =< H andz, «< u,.A 


large uw. implies a small 7. as 





T, « (75) 


and hence a smaller temperature profile gradient. In Equation (75), the overbar refers to 
averaging Over a suitable time period. 

In this experiment, however, the wind speed was relatively low, with average wind 
speed of about 3.95 m/s. The dependence of roughness on wave steepness may be reduced. 


It may suffice to use Charnock’s formula that has been interpolated to cover regime between 
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low wind and high wind speed [Ref. 16], 


u 
PaO oe (76) 
g 


where g is the gravitational acceleration and w. is the friction velocity. The roughness length 


for a smooth surface is given by 


2, = 0.11 — : (77) 


where vis the dynamic viscosity of air. The total roughness length is given by 


Zo, = 2, TF Zo (78) 


The roughness lengths Z,, and z,, are given by 


Vv 
Zor = a (79) 


oa 
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and 


Vv 
20g = ae : (80) 


4 


Secondly, the longer surface waves do not contribute directly to the roughness, z,, but 
they distort the air flow immediately above the ocean surface and are responsible for the 
vertical displacement of the frame of reference. Furthermore, the presence of the longer 
waves causes a redistribution of wave energy and the slope density of the short waves. This 
in turn influences the direction of the mean stress vector by creating an anisotropic roughness 
field at the surface. The long waves strained the shorter waves, causing them to turn in the 
direction of the longer waves. Donelan [Ref. 17] showed that the presence of long, gentle 
swell propagating in the wind direction has , significant effect on the wind sea spectrum. 
There is also evidence that counter swell, propagating against the wind led to higher drag 
coefficients than those experienced in pure wind waves [Ref. 14]. Therefore the presence of 
swell will influence the turbulence of the air above it. By direct interaction with the air flow 
and due to the fact that swells travel faster than the wind, the swells will transfer momentum 
to the air flow above it and thereby enhance turbulent mixing. Since the random motion of the 
surface waves generate turbulence in the air above the surface waves and increase turbulent 
mixing, it is appropriate to employ a statistical mixing length scale that accounts for the 
increased mixing. The solution adopted here is to characterize the increased turbulence by an 


effective random mixing length scale that varies with the significant surface wave height, /7,. 


7G 


The significant wave height is one of the common parameter that oceanographers use as a 
measure for the average height of ocean surface waves. It is defined as the average of the 


highest one-third waves. The vertical profile of € becomes 














EL[ (+H, 
eM] 
and 
oo ee a 
dz K(z+H) 04 :| | ,?) 


Equation (82) indicates that the presence of ocean surface waves will reduce the profile 
gradient in a random manner that is determined by the surface wave height. An increase in 
wave height will lead to a reduction in the profile gradient. 

With the modified temperature profile, it was demonstrated by the ray tracing program 
that optical continuity between the transmitter and the receiver can be established. This is 
shown in Figure 24. This approach is based on applying the presumed scaling for the near 


surface condition. 
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Ray tracing using modified temperature profile, 30 May 96, 1225 hr 
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Figure 24. Ray Tracing with Proposed Temperature Profile. Optical continuity is established. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


This research showed that for optical signal samples that are not occluded by ocean 
surface waves, the signal fluctuations exhibit a log-normal distribution. This is characteristic 
of signals propagating through a random inhomogeneous medium. For signal samples that are 
occasionally obscured by ocean surface waves, the signal fluctuations deviate from a log- 
normal distribution. 

The effects of ocean surface waves on optical signal fluctuations are revealed in the 
temperature fluctuation spectral densities. The temperature spectral densities display 
prominent peaks at the corresponding peak frequencies of the ocean surface waves. These 
characteristic peaks are not present in the temperature spectral densities computed from 
meteorological measurements over land. | 

Optical signal fluctuations increase as the signal ray path gets closer to the ocean 
surface waves demonstrating the significant effects of the ocean surface waves on signal 
fluctuations. Signal fluctuations were also shown to be dependent on wave age and the 
influence of ocean swells. Younger waves are rougher and are likely to cause a greater 
increase in signal fluctuations than more mature waves cause. When heights of signal samples 
are sufficiently low so that the ocean waves occasionally obstruct the ray paths, the ocean 
waves modulate the optical signals. Under moderate or high wind conditions, the ocean wind 
wave steepness should be incorporated into the roughness length scale, z,, . 


On larger scale variations, the ocean surface waves modify the vertical temperature 
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profile and hence the refractive index profile. The standard logarithmic temperature profile 
derived for meteorological conditions over land is not suitable for the marine environment 
above the undulating ocean surface. Ray tracing using the standard logarithmic temperature 
profile indicates that optical continuity could not be established between the transmitter and 
the receiver. Unlike conditions over land, the ocean waves distort the air flow above and 
increase turbulent mixing. This increased mixing is proportional to the ocean wave height. 
The increased turbulent mixing can be accounted for by incorporating a statistical mixing 
length proportional to the ocean wave height. The proposed temperature profile assures 
optical continuity between the transmitter and the receiver. 

The research has demonstrated the significant effects of ocean surface waves on the 
propagation of low altitude optical signals above the ocean. The results of this study provide 
an insight into the effects of ocean surface waves on an optical signal skimming the ocean 
surface from an engineering view point. The study is not intended to disprove existing 
theories and hypotheses but rather to emphasize the relative importance of certain parameters 
that affect the propagation of optical signals in the vicinity of the ocean surface where the 
influence of ocean surface waves is significant. 

The existing model (e.g. IR Tool) used by the Navy utilizes the standard logarithmic 
temperature profile to predict refraction of signals. The standard logarithmic profile is only 
valid for heights that are several significant wave heights above the ocean surface waves 
where the ocean wave effects are negligible. Within a few meters of the ocean surface waves, 


excessive ray bending occurs due to the logarithmic singularity of the temperature profile as 
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the ray approaches the ocean surface. The proposed statistical mixing length approach to 
modifying the temperature profile is perceived as an engineering fix to alleviate the large 
gradient in this critical layer above the ocean surface where sea-skimming threats are 
predominant. Further refinement in the modeling is recommended through extensive 
experimental measurements and data analysis. Modeling of signal propagation above an air- 
ocean interface spans an inter-disciplinary study of electronic engineering, oceanography and 
meteorology. The complexities of the air-ocean interface presents a unique challenge to 
scientists and engineers in a continuous search for more accurate modeling of signal 


propagation in this critical marine layer. 
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APPENDIX A. OPTICAL SIGNAL SAMPLES 


ua, 


Optical Signal Samples 


| | A | |) ee eee eee eee eee 
_1 [Date time |Rxht[m] [Tx ht{m] |height of [rms |% ——[h’2_— | togvvar| _—_ log var 
2t CC frayppath wavet [| | __s wave 60) 
| 3 | 5301225) 4.69] 5.16) 1.71] 0.36] 0.000106) 18.7128] 0.079] 
4 | 5301230) 4.56) 5.18] 1.64] 0.36] 0.000248] 16.5491/ 0.082; 
|S { 5301233} 4.41] 5.19] 1.57 0.36) 0.00062/ 14.9584; 0.095, 
6 | 5301236|  4.27/ ~—S—«5.20/ —=«1.50| ~—=—0«39| (0.005863| 14.5668) 0.067; 
.7 | 5301838] 4.44, 4.94] 1.47/07 1.789256] 10.6063) | (0.229 
| 8 | 5301843; 4.26) 4.91] 1.36) 0.58] 0.956063 7.9751] | (0.198 
| 9 | 5301859; 3.86) 4.82) 1.11) 0.45/ 0.688336] 5.7584; | 0.198 
[4o[ 5301903] 3.77| —4.80| ——1.05| —0.46| 1.123754| 5.1516] «| =—i.145 
(11] 5301909] 3.59] 4.78] 0.94) 0.63] 6.725038] 4.8771] | (0.089 
[12] 5310943) 4.40] 4.95) 1.45] 0.2) 1.78E-11] 44.6558] 0.035, 
113] 5310046; 4.24 4.94] 1.37) 0.19] 3.23E-11) 39.5961] 0.065] 
114] 5310951) 4.15) 4.93] 1.31) 0.22) 1.22E-07| 36.3271] 0.123) 
15] 5311115) 4.17 4.87] 1.30] 0.25] 1.06E-05/ 40.2262) 0.108] 
(16{ 5311120| 4.10) 4.88) 1.27 0.23] 1.88E-06) 37.5196] 0.143] 
47] 5311122) 4.03) 4.88) 1.23] 0.2] 4.1E-08) 35.2933) 0.145] 
18] 5311830] 3.86] 4.87] 1.13) (0.33) 0.030874) 15.0818] 0.154] 
(19| 5311834) 3.69) 4.84) 1.03) 0.32) 0.06726) 11.8792; 0.158] 
20] 5311836) 3.53] 4.83) 0.93) 0.26] 0.016937) 9.2211; 0.157) 
(21| 5311839; 3.43] 4.82) 0.87) 0.34] 0.50895) 8.2193; 0.149) 
(22| 5311932) 2.94 4.55) 0.48) 0.26] 3.26005) 3.2021 | (0.087 
[23| 5311937; 3.07} 4.53) 0.54] 0.29] 3.116286] 3.3635] | (0.294 
[24] 5311948] 3.32) 4.47| 0.66} 0.26) 0.586267) 4.4495] | 0.244 
| 25} 6012208 3.22 4.05 0.41| 0.26) 5.780919| 3.5517/ = | 0.227 
(26| 6012243] 3.21) 4.03] 0.391 0.2] 2.511705) 3.3162; | (0.252, 
| 27| 6012308] 3.53] 4.05] 0.57| 0.19] 0.134146] 6.4756] 0.134) 
(28{ 6021245] 4.391 3.42) 0.87] 0.19 0.010306/ 12.4613) 0.257|_ 


(29 | 6021417 2.43 3.53| -0.23 0.22; 85.2618] 3.0506, __ 0.086 
30| 6021451, 1.84 3.69] -0.52 0.2) 99.52231| 7.2935 [0.103 


/31{ 6030930/ 4.01, 4.01) 0.79) 0.22) 0.015119/ 9.949) 0.107) 
/32{ 6030937; 3.81) 3.96] 0.67; 0.22/ 0.117065) 6.8136) 0.126) 
[33] 6031540) 4.26] 5.19) 1.49) 0.38] 0.004213) 24.3964] 0.103] 
(34] 6031551) 3.83] 5.23] 1.28) 0.3] 0.001032) 17.7704] 0.091, 
[35] 6031605) 3.34) 5.26) 1.02) 0.31| 0.052834) 11.5735] | 0.132 
(36] 6031611; 3.43] 5.28] 1.07) 0.34] 0.079321) 13.1283] || 0.079 
}37| 6041325) 3.58] 4.74) 0.92) 0.26) 0.020864) 11.8525] 0.157) 
(38{ 6041417| 2.77] 4.68) 0.44) 0.34] 9.854486] 3.1197) | 0.269) 
/39{ 6041426] 3.38] 4.69 0.79) 0.28| 0.249417] 7.3583] 0.145) 
|40] 6051933] 4.63] 3.38] 0.76) 0.31] 0.703241] 10.8039) 0.181) 
}41{ 6051945] 3.42] 3.40) 0.20) 0.24) 2057304) 1.5127) | 0.196 
/42| 6052000) 1.69] 3.42) 0.72) 0.19] 99.9921| 8.4996] | 123: 
[43] 6052041) 2.42] 1.69) -1.17/ 0.22) 99.99999] 19.4317] | 0.07, 
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Legend: 


Column A: 


Column B: 


Column C: 


Column D: 


Column E: 


Column F: 


Column G: 


Column H: 


Column I: 


Date and time of signal sample. 

Receiver height above mean sea level (in meters). 

Transmitter height above mean sea level (in meters). 

Lowest height of the ray path above mean sea level (in meters). 

Rms ocean wave height (in meters). 

Percentage of ocean surface waves protruding into optical signal ray path. 
Normalized mean square distance of optical signal ray path above ocean 
surface waves. 

Log-normal variance of signal samples not occluded by ocean surface waves. 
Log-normal variance of signal samples occluded by ocean surface waves. This 
is not strictly correct as the signal samples no longer exhibit log-normal 


distribution. 


i 


Oo Sl U7 


ae 





APPENDIX B. MIXING LENGTH HYPOTHESIS 


The mixing length concept was introduced by Prandtl in 1925 to describe the 
properties of turbulent transfer. In this model, it is assumed that there exists a ieee scale 
called the mixing length over which a turbulent eddy can travel and still maintain its identity. 

Suppose, under neutral atmospheric conditions, a parcel of air moves vertically 
through a small distance, /, from its original position, z+/, where the characteristic quantity, 


q(z), has an average value of g(z+/), then 


q(z) = g(zth) . (83) 


Using Reynold’s decomposition, g(z) can be wnitten as a sum of an average component and 


a fluctuating component, 


Q(z) = gz) + g'() . (84) 


Using mixing length hypothesis, 


q'(z) = q+) - gq) = ed . (85) 


Suppose 2 > 0, then as shown in Figure 25, / > 0 implies that the eddy comes from 
Z 


above and hence the vertical turbulent velocity, w’ < 0. By the same token, / < 0 implies 


8] 


w’ > 0. The above argument indicates that / and w’ are negatively correlated, 1.e., 
(36) 


We 


Assuming that mixing length hypothesis applies to the turbulent velocity component 


and the air temperature then 


(cy ae (87) 
Oz 
and 
(88) 


Cu 
use) az 


where 7’ 1s the turbulent air temperature and wz’ is the turbulent streamline wind velocity. 


Equation (88) implies 
wu’ = w'l (89) 


Wu 


Sas 


In the atmospheric surface layer, < > 0 and since w’/ < 0, then 
Z 


/.,f 


wu <0. (90) 


Assuming the eddy velocity components, uw’ and w’ have the same order of magnitude, then, 


/ ! j0u 
we~ iu’ ~ pen (91) 
Equations (89), (90) and (91) lead to 
— —\ 2 
7! = rf ) | (92) 


In the atmospheric surface layer, the mixing length varies with height above the earth’s 
surface. Near the earth’s surface, mixing cannot take place over large distances due to the 


solid boundary. Away from the surface, the turbulence is less inhibited and mixing is 


enhanced. In general, it is assumed that 


(93) 


where kK is the von Karman constant. 


In the atmospheric surface layer, the main task is to establish relationships between 


83 


vertical gradient of mean quantities of wind speed, temperature and specific humidity and 
surface fluxes. In particular, the surface fluxes of momentum, t, sensible heat, O,,, and latent 


heat, O,, are given by 


t = -pu'w' = pur, (94) 
QO, = pe, Tw’ = -pc,u,T, (95) 

and 
QO, = pL,q'w’ = -pL,u,q, . (96) 


where c, is the specific heat of air, and L, is the latent heat of vaporization of water. 


Combining Equations (92), (93) and (94), 


u 
oe (97) 





where uw is the mean wind velocity (the overbar has been dropped for convenience). Similar 


equations for 7 and g can be written: 
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8 
Oz KZ es) 
and 
oq = ae (99) 
Oz KZ 


For non-neutral atmospheric conditions, Equations (97), (98) and (99) become 


a u, 
3: "Kg PCL) , | one 
or sf, 
ae ,(2/L) (101) 


and 
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og _ 7, 
i ae (Dyn 
oz Kz O,2 ) 


where ’s are the stability functions determined empirically, 


_J(1-202z/L)™, = z/L<0 
1+8z/L, ZL20 ~ 


o- 


b,= 0 = | - 162/L)*?, z/L<O 
‘ ; eZee, z/L>0 


(102) 


(103) 


(104) 


and 2/L <0, 2/L >0, and z/L = 0 correspond to unstable, stable and neutral stratification 


respectively and L is the Obukhov length. 


Integrating equations (100), (101) and (102) yield profiles for mean wu, 7 and q, 


u=u.+—-[In(z/z),)-W,(Z/L) , 
K 
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(105) 


Z | 
T= 1, +— [In (@lzop) -¥7@IL)] (106) 


and 


‘ [In(z/z,,) -p,(2/L)] , (107) 





q=4,* 


where u,, 7, , g are the values at the surface and w’s are the empincally derived stability 


correction functions given in Chapter V. 
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zt+l 


Figure 25. Mixing Length 
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I. TEMPERATURE SPECTRUM AND THE CORRESPONDING WIND WAVE 
SPECTRUM 
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Temperature spectrum, 30 May 96, 1230 hr 
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Ocean surface elevation spectrum, 30 May 96, 1230 hr 
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Temperature spectrum, 30 May 96, 1236 hr 
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Temperature spectrum, 30 May 96, 1838 hr 
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Ocean surface elevation spectrum, 30 May 96, 1838 hr 
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: Temperature spectrum, 30 May 96, 1859 hr 
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Temperature spectrum, 30 May 96, 1909 hr 





Ocean surface elevation spectrum, 30 May 96, 1909 hr 
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Temperature spectrum, 31 May 96, 946 hr 





Ocean surface elevation spectrum, 31 May 96, 946 hr 
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Temperature spectrum, 31 May 96, 1115 hr 
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Temperature spectrum, 31 May 96, 1120 hr 
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Temperature spectrum, 31 May 96, 1830 hr 
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Temperature spectrum, 31 May 96, 1948 hr 
















10 
Sees 
1 OQ a ee 
ee 
i) 
N 
Lt 2 
GN 10 ae EES BE SS ES eee SSS ee) EES se) SS es eS SS GS Ee EE ee 
< ee ee ee eee | lee eee 
©) = 
— 
———— = 
a ee ee tN OT Oe 
SS Se CO a SL + aa ue 
foe ; 4 





= 2 
COC 
10' 


10° 


Ocean surface elevation spectrum, 31 May 96, 1948 hr 


== 
ee alae 


ee 





Temperature spectrum, 1 Jun 96, 2208 hr 
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Ocean surface elevation spectrum, 1 Jun 96, 2308 hr 
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Temperature spectrum, 2 Jun 96, 1245 hr 
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Ocean surface elevation spectrum, 3 Jun 96, 1551 hr 
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Optical signal spectrum, 2 Jun 96, 1451 hr 
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Ocean surface elevation spectrum, 5 Jun 96, 1933 hr 
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IV. RAY TRACING USING STANDARD LOGARITHMIC TEMPERATURE 
PROFILE 
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Ray tracing using standard logarithmic temperature profile, 30 May 96, 1225 hr 
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Ray tracing using standard logarithmic temperature profile, 30 May 96, 1230 hr 
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Ray tracing using standard logarithmic temperature profile, 30 May 96, 1233 hr 
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Ray tracing using standard logarithmic temperature profile, 30 May 96, 1236 hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 943 hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 946 hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 951 hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 1115hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 1120 hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 1122 hr 


Source ht = 4.88 m 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 1830 hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 1834 hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 1836 hr 
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Ray tracing using standard logarithmic temperature profile, 31 May 96, 1839 hr 
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Ray tracing using standard logarithmic temperature profile, 1 Jun 96, 2308 hr 
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Ray tracing using standard logarithmic temperature profile, 3 Jun 96, 930 hr 
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Ray tracing using standard logarithmic temperature profile, 3 Jun 96, 937 hr 
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Ray tracing using standard logarithmic temperature profile, 4 Jun 96, 1325 hr 


Source ht = 4.74 m 
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Ray tracing using standard logarithmic temperature profile, 4 Jun 96, 1426 hr 


Source ht = 4.69 m 


& 
CO 
O £ 
© <x 
WW = © 
a 
tL 

@ iil 
2 
o oD 
oS «= 
oOo Ow 
Yr wx 


= 
c€ 
pe 
© 
Uf 
~~ 
i= 
® 
> 
S 





uw wT 
[i] yUBIS}H 


Range [km] 


V. RAY TRACING USING MODIFIED TEMPERATURE PROFILE 
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Ray tracing using modified temperature profile, 30 May 96, 1225 hr 
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Ray tracing using modified temperature profile, 30 May 96, 
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Ray tracing using modified temperature profile, 31 May 96, 946 hr 
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Ray tracing using modified temperature profile, 31 May 96, 951 hr 
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Ray tracing using modified temperature profile, 31 May 96, 1115 hr 
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Ray tracing using modified temperature profile, 31 May 96, 1120 hr 
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Ray tracing using modified temperature profile, 31 May 96, 1122 hr 
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Ray tracing using modified temperature profile, 31 May 96, 1830 hr 
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Ray tracing using modified temperature profile, 31 May 96, 1834 hr 
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Ray tracing using modified temperature profile, 31 May 96, 1836 hr 
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Ray tracing using modified temperature profile, 31 May 96, 1839 hr 
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Ray tracing using modified temperature profile, 1 Jun 96, 2308 hr 
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Ray tracing using modified temperature profile, 3 Jun 96, 930 hr 
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Ray tracing using modified temperature profile, 3 Jun 96, 937 hr 
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Ray tracing using modified temperature profile, 4 Jun 96, 1426 hr 
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